Bose-Einstein Condensates in Optical Lattices: Band-Gap Structure and Solitons 



<N : 
o . 
o ■ 
(N ; 

3 ■ 

<: 

t> ■ 

(N ■ 



> 

in 
oo 

o 

CM 
O 



C 

O 

o 



13 



Pearl J.Y. Louis 1 , Elena A. Ostrovskaya 1,2 , Craig M. Savage 2 , and Yuri S. Kivshar 1 
1 Nonlinear Physics Group, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 
2 Department of Physics and Theoretical Physics, 
Australian National University, ACT 0200, Australia 

We analyze the existence and stability of spatially extended (Bloch-type) and localized states of 
a Bose-Einstein condensate loaded into an optical lattice. In the framework of the Gross-Pitaevskii 
equation with a periodic potential, we study the band-gap structure of the matter-wave spectrum in 
both the linear and nonlinear regimes. We demonstrate the existence of families of spatially localized 
matter-wave gap solitons, and analyze their stability in different band gaps, for both repulsive and 
attractive atomic interactions. 
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I. INTRODUCTION 

Trapping of neutral atoms in optical lattices Q| pro- 
vides an effective and powerful means for controlling 
the properties of Bose-Einstein condensates (BECs), as 
confirmed by a number of the recent experiments (see, 
e.g., Refs. jg § | |, |, % |). Coupling between the 
BEC droplets localized in weakly interacting minima of a 
standing light wave creates a reconfigurable matter-wave 
structure that can be easily manipulated by varying the 
strength of the periodic potential and its well spacing. 

Many properties of such arrays of the BEC droplets 
can be understood, within the framework of the mean- 
field approach, by employing the tight-binding approx- 
imation borrowed from solid state physics (|, [To[ pd| . 
Being based on the assumption of strong localization of 
the condensate wave functions in the individual potential 
wells of the lattice, the tight-binding approximation and 
the resulting discrete lattice models provide a limited de- 
scription of the BEC dynamics in an optical lattice. An 
alternative reduction of the continuous system to the dis- 
crete vector lattice involves employing a complete set of 
on-site Wannier states |l2], and is more general than a 
tight-binding approximation. However, it is to be clari- 
fied whether such a reduction can provide predictions on 
the dynamics and stability of a BEC in a periodic po- 
tential more readily than the full mean-field model. A 
more adequate analysis of both periodic and localized 
modes of the continuous Gross-Pitaevskii (CP) equation 
with a periodic lattice potential is still missing. Some of 
the recent studies of the GP model [Hj, p|, [E§] point 
towards interesting features of the BEC dynamics that 
should be observed beyond the applicability of discrete 
models. 

The main purpose of this paper is to analyze the struc- 
ture and stability of the stationary states of a condensate 
in an optical lattice, the so-called nonlinear modes |fl7| , 
in the framework of a continuous mean-held model. This 
approach allows us to consider more general types of sta- 
tionary states and BEC dynamics without many of the 
restrictions of the discrete models. We solve the contin- 



uous one-dimensional GP equation numerically and de- 
scribe families of stationary states, which possess differ- 
ent spatial structure and stability properties. 

First, we find spatially extended stationary nonlinear 
periodic modes of the condensate in a lattice. These pe- 
riodic modes possess different symmetry, and provide a 
generalization of both the Bloch modes in the theory of 
linear periodic systems and of the extended states of a 
condensate in a finite optical lattice with a confining po- 
tential jl8| . We demonstrate that the matter- wave spec- 
trum for such periodic BEC states has a band-gap struc- 
ture, which survives in the presence of the mean-field 
nonlinearity in the GP equation. 

We also find novel types of condensate states which 
cannot exist in the linear limit. These are spatially local- 
ized modes of the condensate which exist in the band gaps 
of the matter-wave spectrum and, therefore, they should 
be called gap solitons. The existence of the matter-wave 
gap solitons was first suggested in the framework of the 
coupled- mode theory Q , in analogy to optical gap soli- 
tons in Bragg gratings. Remarkably, gap solitons - spa- 
tially localized atomic wavepackets - may exist even in 
a condensate with repulsive interatomic interactions p3[ . 
Their experimental demonstration, however, remains a 
challenge. Here, we present a comprehensive theory of 
the structure and linear stability analysis of different 
families of gap solitons with respect to variations in the 
lattice parameters and effective nonlinearity defined by 
the number of atoms in the condensate. We discuss the 
key scenarios of instability of spatially localized modes in 
optical lattices and the physical mechanism for creating 
different multiple states of the BEC droplets. 

We would like to emphasize that many of the effects 
described here in the framework of the GP equation with 
a periodic potential, such as the mode stability and the 
existence of the localized modes in different band gaps, 
cannot be obtained in the framework of the tight-binding 
approximation since the coherent and spectral properties 
of these two models arc different. This latter issue is simi- 
lar to the BEC dynamics in a double- well potential where 
a rigorous coupled-mode wave theory should be based on 
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the analysis of nonlinear modes of the entire structure, 
rather than condensate wave functions of isolated wells 
& 

The paper is organized as follows. In Sec. II we present 
our model and reduce it to the one-dimensional GP equa- 
tion, under the assumption of a cigar-shaped conden- 
sate. Section III includes the analysis of both linear 
and nonlinear matter- wave spectra in an optical lattice. 
We revisit the structure of the band-gap spectrum for 
the Bloch-like modes of the noninteracting condensate, 
and extend this concept to the nonlinear case, demon- 
strating the existence of the nonlinear Bloch-like modes 
and the nonlincarity-dependent band-gap structure of the 
matter-wave spectrum in optical lattices. The study of 
the BEC gap solitons in different band gaps is summa- 
rized in Sec. IV, where we also present our results on 
the stability of gap solitons, for both repulsive and at- 
tractive atomic interactions. Section V concludes with a 
summary and further perspectives. 



II. MODEL 

We consider the macroscopic dynamics of a BEC 
loaded onto an elongated optical lattice, similar to those 
created in a number of experiments (see, e.g., Refs. 
jl], ||, ||). The condensate dynamics is described by the 
three-dimensional GP equation, 
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where ^(X,Y, Z;T) is the wave function of the cigar- 
shaped condensate, (Y, Z) are the directions of strong 
transverse confinement, and X is the direction of the 
lattice. The combined potential of the optical lat- 
tice and magnetic trap V(X, Y, Z) can be written as, 
V{X,R 2 ) = E sm 2 {TrX/d) + (l/2)m(uj x 2 X 2 + to\R 2 ), 
where R 2 = Y 2 + Z 2 , E is the well depth of the opti- 
cal lattice, d is the characteristic lattice constant, and 
are trapping frequencies in the corresponding directions. 
The parameter T = 4irh 2 a s /m characterizes the s-wave 
scattering of atoms in the condensate which introduces 
an effective nonlinearity in the mean-field equation; it is 
positive for repulsive interactions and negative for attrac- 
tive interactions. 

Measuring time in units of wj 1 , the spatial variables 
in units of the transverse harmonic oscillator length, 
ao = (h/mLu^) 1 / 2 , the wave function amplitude in units 
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and the potential in units of hujj_ , we obtain the 



following dimcnsionless GP equation: 
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where t = Tuj±, (x,r) — (X,R)/ao, and <?3D = 
47r(a s /do) . The potential now takes the form V(x, r 2 ) = 
V sm 2 (Kx) + (l/2)(tt 2 x 2 + r 2 ),where V = Eq/{Twj x ), 
and K = ttoq / d. The ratio of the confinement strengths 



for the magnetic trap £1 = ui x /uj± varies from 10 1 
§,|,|,[§tol/V2§0§. 

To simplify our analysis, we consider a cigar-shaped 
condensate in a strongly elongated trap (f2 ~ 1CP 1 ) 
and a quasi-one-dimensional optical lattice in the di- 
rection of a weak confinement. Under these assump- 
tions, the model can be reduced to a one-dimensional GP 
equation by assuming separable solutions of the form, 
^(r, x;t) = $>(r)ip(x;t). Due to the tight confinement 
of the condensate in the transverse direction, we as- 
sume that the spatial structure of the function <f>(r) is 
well described by a solution of the two-dimensional ra- 
dially symmetric quantum harmonic oscillator problem, 
Vj_3> + $ — r 2 $ = 0. Taking the ground state solution, 
$(r) = Ccxp(— r 2 /2), where C 2 = 1/n is found from 
the normalization condition, /_ \$>\ 2 dydz = 1, and inte- 
grating Eq. ([!]) over the transverse coordinate, we derive 
an effective one-dimensional GP equation in the dimen- 
sionless variables, 
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where (?id = 2\a 8 \/ao, and the coefficient a — sgn(a s ) = 
±1 characterizes the type of the s-wave interactions. Due 
to the weak magnetic confinement in x direction, it is 
possible to omit the term f2 2 from the effective po- 
tential, thus assuming that the ID lattice potential , 
V[x) — Vo sin 2 (i^cc), is uniform and is not affected by 
the presence of the additional confining potential. Fur- 
thermore, it is convenient to reduce the parameter space 
of the model i.e. effectively set <?id = 1) by rescaling 
the wave function as ip — * ijj^/giu- In physical terms, it 
means that the peak density of the condensate is used to 
characterize the effective mean-field nonlinearity. 

The model (^) can describe the BEC in a one- 
dimensional optical lattice over a wide range of exper- 
imental parameters, such as lattice depth, well spacing, 
and the effective nonlinearity. Indeed, Eq. (||) is invari- 
ant with respect to the following scaling transformation: 
{K,V ,x,t,^} -> {K/i],Vo/ri 2 7 x7j,tr] 2 ,il}/r]}, where rj is 
a dimensionless constant. However, a significant point 
of difference between the lattice in our model and those 
in experiments jl], 0, |3| is the absence of an additional 
confining potential (e.g., gravity or magnetic trapping). 

We determine the experimentally relevant range of di- 
mensionless parameters for our model by taking the case 
of 87 Rb atoms, which sets m — 1.44 x 1CT 25 kg and 
a s = 5.3 nm. The characteristic lattice constant, d, is 
determined by the angle between the intersecting laser 
beams forming the lattice and, in the current experi- 
ments, can be varied in the range 0.4 — 1.6 /im || 
The lattice depth, Eo, scales linearly with the light inten- 
sity, and varies between zero and £'™ ax ~ 20E 1CC , where 
E mc = h 2 n 2 /(2md 2 ) is the lattice recoil energy 0, 
Taking the transverse frequency as luj_ ~ 2tt x lCr Hz 
H Q , we obtain the following regions of the model param- 
eters: 2 < K ^ 8 and 4.5 x 10 1 V£ lax ^ 7.0 x 10 2 , re- 
spectively. In what follows, for the sake of computational 



3 



convenience, we will make use of the scaling properties 
of the model and use two different sets of the parame- 
ters: K = 0.1, V ~ 0.1 (ry = 5), and K = 0.4, V ~ 1.0 
(n = 20). 

We note that the number of atoms in the condensate 
wavefunction, TV, is now given by the expression M = 
N/giD, where 



N = 



\ip(x,t)\ 2 dx 



(3) 



is an integral of motion for the model 



III. BAND-GAP SPECTRUM OF MATTER 
WAVES 

A. Linear Bloch modes 

Stationary states of the condensate in an infinite peri- 
odic potential of an optical lattice ('nonlinear modes') 
are described by solutions of Eq. (||) of the form: 
ij)(x,t) = 4>{x) exp(— ifit), where \i is the corresponding 
chemical potential. The steady-state wave function 4>{x) 
obeys the time-independent GP equation: 

^S" [^sin 2 (X*)0-^] -<tM 2 = O, (4) 

where the lattice depth Vq and period K are the exper- 
imentally controlled parameters, and they both can be 
varied in our model. 

In the case of noninteracting condensate, which for- 
mally corresponds to a — 0, Eq. (Q) is linear in tf> and, 
as has been shown in [ 21|, the Bloch- function formalism 
applies. In this case (0) can be written as the Mathieu's 
equation: 



i 

dy 



- z |-[2gcos(2»)-p]0 = O, 



(5) 



where y = Kx, q = -V /(2K 2 ), and p = 2(q + ll/K 2 ). 
By employing the Floquet-Bloch theorem, the conden- 
sate wave function can be presented as a superposition 
of Bloch waves, 



4>(y) = biPi(y)e iXy +b 2 P 2 (y)e 



-i\y 



(6) 



where P\^{y) are periodic functions with the least period 
7r, b\ t 2 are constants, and A is the characteristic of the 
Floquet exponent. The structure of all solutions to the 
Mathieu's equation can be determined by Floquet theory 
(this text-book analysis can be found, e.g., in Ref. p2j). 
According to this theory, the spectrum of Eq. (||) consists 
of domains ('bands') of the (p, q) space in which there 
exist only amplitude-bounded, oscillatory solutions. The 
bands are separated by the regions ('gaps') where un- 
bounded solutions exist. The band edges correspond to 
exactly periodic solutions. 
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FIG. 1: Band-gap diagram for Bloch waves in the linear 
regime (noninteracting BEC in an optical lattice). Shaded 
areas: bands of bounded oscillatory solutions for the matter 
waves; numbered areas (I and II)- the lowest two band gaps 
in the spectrum where some solutions are unbounded and no 
periodic solutions exist. Solid lines - band edges correspond- 
ing to periodic solutions. Dashed line shows the particular 
value Vb = 0.1 used below. 



Figure |l| presents the band-gap diagram for the ex- 
tended solutions of Eq. (||) which describe noninteract- 
ing condensed atoms in an optical lattice. The results 
are presented for the parameter domain (/i, Vb) relevant 
to our problem. The shaded areas correspond to the re- 
gions where the Floquet exponent has a real part (which 
is an odd or even integer) and, therefore, spatially oscil- 
lating (but not ingeneral periodic) solutions to the linear 
model equation ffl) exist. If A has a nonvanishing imag- 
inary part, the solutions do not oscillate (decay or grow 
exponentially with x); such solutions belong to the gaps 
separating the bands of oscillatory solutions. The first 
two gaps where unbounded solutions exist, numbered I 
and II in Fig. [I] are bordered by the band edges repre- 
senting the linear Bloch waves with periodicity of 2n/K 
and tt/K, respectively. 

It is in these gaps of the linear spectrum that the non- 
linear localization of matter waves in the form of gap 
solitons occurs, and it is this localization that could en- 
able the creation of solitons in the BEC with repulsive 
interatomic interaction. The chemical potentials corre- 
sponding to such localized waves lie in the band gaps of 
the spectrum of the noninteracting condensate (marked 
areas in Fig. Q). 

The dashed line in Fig. [I] marks the value of the lattice 
depth for which the the study of the gap modes presented 
below is conducted. We note that this value is greater 
than the threshold value V* — fx that roughly marks the 
transition between the quasi- unbounded (Vb < V*) and 
strongly bounded (Vb > V*) condensate wavefunctions 
(see also discussion in Ref. p3fl). The tight-binding or 
"fragmented" condensate regime corresponds to Vb ^> 
V* , where the bands " collapse" to discrete levels of the 
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Wannier states pjj . 

The knowledge of a complete band-gap spectrum of 
the matter-waves in an optical lattice provides us with 
important clues on the existence and stability regions 
of different types of BEC localized states that may be 
excited in an optical lattice. The frequently employed 
nearest-neighbor tight-binding approximation based on 
a discrete model jlOl [ll], 24 is inferior to the analysis 
of the complete continuous model in that it describes 
only one (the first) band surrounded by two semi-infinite 
gaps. The following sections aim to provide us with un- 
derstanding of the details of the formation and stability 
of "gap modes" - the nonlinear localized states of the 
condensate in a lattice. Despite the recent advances in 
the treatment of the continuous GP models |l4|, [ll], ||| , 
such an analysis is still missing |38[. 
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FIG. 2: Middle: Condensate population (shown as the num- 
ber of atoms, N, denned by Eq. ([|), per lattice period, tt/K) 
for the extended stationary states of a repulsive interaction 
(a > 0) BEC in an optical lattice at Vo = 0.1. The values 
of n corresponding to the linear band edges at Vo = 0.1 are 
recovered in the limit N — » (cf. Fig. h dashed line). Top 
and bottom: dotted lines - structure of the lattice potential, 
and solid lines - spatial structure of the nonlinear Bloch waves 
at the marked points on the band edges. 



B. Nonlinear Bloch modes 

The effect of interatomic interaction, described by the 
cubic nonlinearity in the mean-field GP equation, results 
in an effective shift of the chemical potential for the ex- 
tended periodic states, i.e. Bloch waves of BEC, on the 
band edges. This leads to the corresponding shifts of the 
band-gap edges, and the gaps can even close up for large 
densities of the repulsive condensates. The nonlinearity- 
induced shift of the bands can be described by the multi- 
scale perturbation theory developed for the Mathieu's 
equation with cubic nonlinearity (H) and a shallow po- 
tential (Vo <C 1) (see, e.g. Ref. |p6| ). Here, we study 
the extended periodic stationary states of the interact- 
ing condensates in the form of nonlinear Bloch waves, 
by solving Eq. (Q) numerically with standard relaxation 
techniques. 

In correspondence to the linear odd and even Bloch 
waves that satisfy Eq. (||), the nonlinear Bloch waves 
display similar spatial structures, in which the droplets 
of the neighboring wells are out of phase or in phase with 
each other (see Fig. ||). The families corresponding to 
the two distinct stationary states for different condensate 
densities are shown in Fig. ||for the repulsive condensate 
((7 = +1). The departure of the nonlinear Bloch states 
from the linear limit (recovered at N — > 0) indicates the 
degree to which the mean-field nonlinearity affects the 
structure of the band gaps in the matter-wave spectrum. 
However large this departure is, the spatial periodicity 
and symmetry of the linear Bloch waves corresponding 
to different bands (see, e.g., Ref. pTj) is preserved in 
the nonlinear regime. The " in-phase" and " out-of-phase" 
extended periodic states of the condensate corresponding 
to the lower and upper edges of the lowest band (see Fig. 
||, bottom panel) represent a continuous-wave analogue 
to the so-called " unstaggered" and " staggered" modes of 
the discrete lattice model, respectively However, 
none of the higher-order band modes can be described 
by the discrete model, in the framework of the nearest- 
neighbor tight-binding approximation, and its validity is 
therefore restricted to the lowest band only. 



IV. MATTER- WAVE GAP SOLITONS 

A. Spatially localized modes 

We have numerically found different families of spa- 
tially localized BEC states in an optical lattice, for both 
repulsive and attractive interaction. Such states can be 
regarded as strongly localized soliton-like modes created 
in band gaps of the optical-lattice spectrum, and there- 
fore they should be compared with gap solitons found for 
other periodic systems beyond the coupled-mode theory 
(see, e.g., Refs. [^8[ |9| and references therein). Several 
families of the lowest-order gap modes and their spatial 
structure are presented in Figs. ||| to [?] for a = +1. The 
gap solitons exist in all band gaps, excluding the semi- 
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infinite gap of the spectrum below the first band, which 
is analogous to the total internal reflection gap in optical 
Bragg gratings. 
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FIG. 3: Families of different gap solitons. The structure of 
the modes corresponding to the marked points is shown in 
Figs. ^| , ^, and |(j Shaded areas show the linear Bloch-wave 
bands. 
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FIG. 4: Spatial structure of the localized states in the first 
band gap corresponding to the marked families A, B, and C 
in Fig. 

One of our main results is that each of the band gaps 
includes a branch of the lower-order fundamental gap 



FIG. 5: Spatial structure of the localized states in the first 
band gap corresponding to the marked families D, E, F, and 
G in Fig. 



modes possessing a single central peak. These funda- 
mental localized states, or the BEG droplets, are strongly 
localized near a minimum of a potential well, as shown 
for the first two gaps by the branches A and H in Fig. 
H and the corresponding profiles in Figs. ||, ||, and ||. 
Higher-order modes can be thought of as bound states 
of several fundamental modes, including the continuous 
counterparts of "twisted modes" of the discrete model 
p9| and multiple states of the localized modes in nonlin- 
ear lattices pof . 

In the first gap (see Figs. ^ and || ), we show only 
the lowest-order bound states of two (states B, D, F, 
and G) or three (states C and E) in-phase (states D, E, 
and F) or out-of-phase (states B, C, and G) fundamen- 
tal modes. These can be "even", i.e. centered between 
potential wells (states B, C, D, and E) or "odd" , i.e. cen- 
tered on the potential minima (states F and G). For low 
condensate densities, the number of atoms in the local- 
ized states that correspond to the bound states of two 
or three fundamental modes are, correspondingly, two or 
three times larger than that in a fundamental state (see 
Fig. ^) . The families of the higher-order formed of the 
same number of fundamental states are degenerate in N 
for low densities; this degeneracy lifts for larger densities, 
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FIG. 7: Spatial structure of the states in the second trans- 
mission band, localized on the extended background and cor- 
responding to the families (a) A and (b) B in Fig ^j. 




FIG. 8: Top: the family diagram and bottom: spatial struc- 
ture of the lowest-order localized states of the attractive 
BEC (a — —1) in the semi-infinite gap: (a) far from, and 
(b) near the edge of the first band, for /j, = 5 x 10~ 4 and 
H — 2.22 x 10~ 2 , respectively. 



FIG. 6: Spatial structure of the localized states in the second 
band gap corresponding to the marked families in Fig. ^. 



i.e. for larger mean-field nonlinearities. It is interesting 
to note that, along with the fundamental and its bound 
states, there exist higher-order localized modes and their 
corresponding bound states that do not have analogues 
in the discrete lattice models. A representative mode (K 
and its out-of-phase bound state (L) are shown in Fig 

The spatially localized modes exist only in the gaps of 
the Bloch-wave spectrum. Closer to the band edges, the 
mode structure changes dramatically. In particular, all 
the localized modes of BEC with repulsive interaction ( 
a = +1) develop extended "staggered" oscillating tails 
near the upper gap edges, as demonstrated for the fun- 
damental mode and two of its bound states in the second 
gap (Fig. |[ right column at fi = 0.09). The spatial 
symmetry and periodicity of these tails corresponds to 
the Bloch-wave structure at the corresponding band edge 
(see Fig. |6| second column). Inside the bands, these lo- 



calized states (two examples are shown in Fig. Q) exist as 
the modes on an extended, oscillating background with 
the structure of the corresponding Bloch wave [cf. Figs, 
fj] and H(c)]. These states, localized on a background, 
are analogous to the "antidark" optical gap solitons de- 
scribed in Ref. ]3l| , and can, in principle, be observed 
experimentally. 

For the attractive interatomic interaction (a = — 1), 
additionally to the gap modes described above, novel 
types of localized modes and their bound states are found 
in the semi-infinite gap, i.e. below the lowest edge of the 
first band of the matter-wave spectrum (see Figs. [I] and 
||). The modes of this type, shown in Fig. || are sup- 
ported by the attractive interaction, and they are similar 
to conventional bright solitons with exponentially decay- 
ing tails that have no nodes and are weakly modulated 
by the lattice potential. For any fixed value of the lattice 
amplitude, Vq, with growing chemical potential, the lo- 
calized modes approach the lower edge of the first band 
of the linear spectrum and develop oscillating tails with 
the structure defined by the structure of the lowest-order 
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B loch waves [cf. Fig. 0(a)]. 



B. Stability of localized modes 

Stability of extended and localized modes in nonlinear 
systems is a very important issue, since only dynamically 
stable modes are likely to be generated and observed in 
experiments. To determine the stability properties of the 
localized modes, we consider small perturbations to a 
solution of the GP equation (||) in the form: 



ip(x, t) = tfi(x)e 



-ifit 



+ e 



u{x)e ll3t + w*(x)e 



-i[3*t 



-ifit 



(7) 

where £<1, and <fi(x) is the steady state localized non- 
linear mode. We linearize Eq. (|J) around the localized 
solution and obtain, to the first order in e, the linear 
eigenvalue problem for the perturbation modes: 



(7 



C -o4> 2 
*2 r 

—l-Q 



(3 



(8) 



where 



Co 



2dx 2 



-V{x)-2a\4>\ 2 + /!. 



In this presentation, the modes describing the develop- 
ment of instability have either purely imaginary or com- 
plex eigenvalues (3; in this latter case the instability is 
called an oscillatory instability. We solve the eigenvalue 
problem (||) for the condensate localized states found 
numerically in the previous section by standard finite- 
difference methods or as a matrix eigenvalue problem 
with Fourier or polynomial interpolants for the differ- 
ential operators. 
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FIG. 9: A typical spectrum of the linear eigenvalue problem 
for the higher-order mode of the repulsive BEC (Vo = 1.0, 
K = 0.4), below (left) and above (right) the instability thresh- 
old. Shaded areas indicate linear Bloch bands. 

It is well known from the theory of optical gap solitons 
in nonlinear Bragg gratings, that the linear stability anal- 
ysis of gap solitons is not trivial and presents significant 
computational difficulties |32L |33|, |34|, |35[ . We have con- 
firmed that the same difficulties arise in the stability anal- 
ysis of matter-wave gap solitons. The reason for this is 




FIG. 10: Top: Example of the resonance between one of the 
perturbation modes and the Bloch-wave band edge leading 
to an oscillatory instability of the odd localized mode of the 
repulsive BEC (Vo = 1.0, K = 0.4). Bottom: (a) mode profile 
and (b) its temporal evolution above the instability threshold 
(fj, = 0.65, fi ± Re/? = 0.733, Im/3 = ±0.019). 



the occurence of weak oscillatory instabilities of gap soli- 
tons, which are associated with poorly localized pertur- 
bation eigenfunctions. Detection of such instabilities for 
a finite continuous system, whether by finite-differences 
or pseudo-spectral methods, is highly sensitive to the 
boundary conditions, the size of the spatial domain, and 
thus to the periodicity of the system. In addition, the dis- 
cretization of the continuous system of the spatial extent 
L with M grid points has been shown to produce spurious 
unstable eigenvalues. In general, the convergence of the 
perturbation eigenvalues as {L.M/L} —* oo is exponen- 
tially slow |32| [53). The effective and accurate detection 
of oscillatory instabilities is possible for L — > oo, and the 
infinite domain boundary conditions can be simulated, 
for example, by expanding the wavefunctions in the ba- 
sis of Bloch waves (or generalized plane waves in optical 
case [[m], [35|). The accuracy of this method is limited by 
the number of Bloch functions in the expansion. Alter- 
natively, the Evans function method used, e.g., in pi[ |2^] 
is capable of treating exact boundary conditions at the 
edges of the spatial domain. 

One should bear in mind, however, that the stabil- 
ity results obtained for the infinite domain may not be 
confirmed by the dynamical simulations of the time- 
dependent GP equation due to their essentially finite spa- 
tial domains. (For instance, optical gap solitons found 
to be stable on an infinite domain can be unstable on 
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a finite domain fl33|).The problem is not only a numeri- 
cal one since the physical extent of the optical lattices is 
always limited by the size of the external trapping poten- 
tial. Below, we present stability analysis for a large but 
finite lattice without an external potential. The virtue of 
this analysis is in its qualitative predictions with regard 
to the types and scenarios of the instability that is exhib- 
ited by the localized modes of the BEC in a ID lattice. 
Quantitatively accurate, experimentally relevant calcula- 
tions of both localized states and their stability should 
ideally treat a finite lattice with the parabolic trapping 
potential. 

A typical structure of the eigenvalue spectrum is pre- 
sented in Fig. [9| Since the perturbation eigenfunctions 
are supported by the periodic potential of the lattice, 
the full spectrum of the eigenvalue problem (||), i.e. the 
values of /i ± /3, has the band-gap structure identical to 
that presented in Fig. [I]. In addition, for any fixed /z, 
there exist modes localized inside the gaps of the con- 
tinuous spectrum, the so-called " internal modes" . These 
linear gap modes exist due to the distributed "defect" 
in the lattice potential introduced by the nonlinear lo- 
calized states. As the chemical potential grows, the in- 
ternal modes of the localized states move towards the 
upper edges of the continuum bands (see Fig. ^|) and 
generate an oscillatory instability through the resonance 
mechanism previously described in a different physical 
context pg[ |. Because different "internal" perturbation 
modes with different symmetries can cross into the band 
and resonate with the linear Bloch waves, the instability 
can trigger various types of spatial dynamics, including 
the symmetry breaking instability. 

Applying the stability analysis to the modes identi- 
fied in Sec. IV A, we found that all of the higher-order 
gap modes identified for our model exhibit oscillatory in- 
stability in the regions of the existence domain near the 
upper band edges. For the parameters used in our calcu- 
lations the strong exponential instability associated with 
a purely imaginary /3 does not occur has not been de- 
tected for gap modes. A typical evolution of the internal 
perturbation mode in a gap and birth of the oscillatory 
instability for large /i is shown in Fig. [l(] (top), for the 
even, out-of-phase gap mode. The time scale for devel- 
opment of such weak instabilities is large, which means 
that oscillatory unstable modes can be dynamically sta- 
ble on the scales comparable with the typical lifetime of 
the condensate. 

For the lowest-order fundamental gap modes, our anal- 
ysis has shown the existence of internal modes which, if 
excited, may be responsible for long-lived oscillations. It 
also produced very narrow regions of extremely weak os- 
cillatory instability (Im/3 ~ 10~ 4 — 10~ 5 ) in the close 
vicinity of the band edges. The accuracy of these calcu- 
lations is low since the gap mode is poorly localized in 
this region. As the computed values Im/3 — ■» with in- 
creased accuracy of the method, it is most likely that this 
oscillatory instability is unphysical and the fundamental 
gap solitons are linearly stable. 



In contrast, some of the higher-order localized modes 
of the attractive condensate (a = — 1), in the semi- infinite 
gap below the first band, have been found to experi- 
ence the strong (exponentially growing) linear instabil- 
ity, characterized by a pair of purely imaginary eigen- 
values of the perturbation modes (3. These modes origi- 
nate from neutral modes of the excitation spectrum (with 
Re/3 = 0). An example of such instability development 
via bifurcation of the neutral mode into the imaginary 
plane is presented in Fig. [ll] for an odd bound state of 
fundamental bright solitons. Remarkably, the instabil- 
ity growth rate drops as the mode moves away from the 
lowest band edge. This is due to the fact that the con- 
densate density grows and, ultimately, strong localization 
of large-amplitude BEC droplets occurs at the individual 
wells. In this limit the lattice effects are weak, and the 
property of a localized state is described by the free-space 
GP equation, for which the (inverted for our choice of fi) 
Vakhitov-Kolokolov stability criterion states that the fun- 
damental bright solitons of attractive BEC are stable as 
long as —dN/dfj, > 0, i.e. in the entire semi-infinite gap 
(see Fig. §, top). 




0.2 



FIG. 11: Top: Growth rate of the strong linear instability 
for an even out-of-phase localized mode of the attractive BEC 
(Vb = 1.0, K — 0.4, a — —1); dotted and solid lines corre- 
spond to the double-peaked even in- and out-of-phase pertur- 
bation modes, respectively. Bottom: (a) the mode profile and 
(b) its temporal evolution; the mode is initially perturbed by 
off-setting the position of its center (/i = 0.1, f3 = 0.157i). 

In order to check the results of the stability anal- 
ysis of the stationary solutions for the BEC modes, 
we have employed direct numerical simulations of the 
time-dependent equation using a split-step Fourier 
method which preserves the normalization of the con- 
densate wave function with a high accuracy. To confirm 
the stable dynamics of gap solitons, the results were ob- 
tained by running the simulations for 8000 dimensionless 
time units which is equivalent to a time of 10.6 seconds. 
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Barrett et al. |3(J who created the first BEC directly in 
an optical trap, found that for a condensate of 3.5 x 10 4 
atoms the 1/e lifetime was about 3.5 seconds. However, 
Stamper-Kurn et al. [ |37| observed that in a purely opti- 
cal trap the 1/e lifetime at low densities was longer than 
10 s, although at that timescale significant loss occurs as 
the population goes from 10 7 atoms to 10 5 atoms. The 
results of our dynamical simulations showing the differ- 
ent instability developments for oscillatory and linearly 
unstable higher-order states of repulsive and attractive 
condensates are presented in Figs. [l0](b) and [n](b), re- 
spectively. 

V. CONCLUSIONS 

We have analyzed the stationary states and stability of 
the condensate in an optical lattice, in the framework of 
the GP equation with a periodic potential. We have con- 
firmed some of the earlier predictions for the condensate 
structure based on the effective tight-binding model, as 
well as revealed a number of novel features not described 
by the effective discrete models. In particular, we have 
extended the concept of the Bloch waves to the case of 
a nonlinear model and described the families of the non- 
linear periodic states which generalize the corresponding 
Bloch-wave modes of linear periodic systems. We have 
demonstrated that the spectrum of such periodic modes 
possesses a band-gap structure, which depends on the 
number of atoms, i.e. the effective nonlincarity in the 
mean-field theory. 

In the gaps of the matter-wave spectrum, we have nu- 
merically found the families of novel spatially localized 



modes of the condensate in an optical lattice trapped 
by a few lattice minima. Such modes should be com- 
pared with optical gap solitons, and they can exist even 
for the condensates with repulsive interaction. We have 
analyzed the structure and, more importantly, the lin- 
ear stability of the gap solitons which belong to different 
band gaps, for both repulsive and attractive interatomic 
interactions. We believe that these novel types of spa- 
tially localized, coherent excitations of the condensate 
can potentially be observed in experiments. 

Our studies extend the analysis of the BEC in an op- 
tical lattice beyond the applicability limits of the fre- 
quently used discrete models based on the tight-binding 
approximation, as well as they call for the future exten- 
sive analysis of the condensate stability in optical lattices 
described by the three-dimensional GP equation with a 
periodic potential. In particular, similar to the conden- 
sates in a parabolic trap, we expect that optical lattices 
can support dark solitons, excited on the different types 
of the nonlinear Bloch waves. The stability of such dark 
solitons should differ dramatically from the stability of 
the conventional dark solitons. These studies are now in 
progress. 
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